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Abstract 

The electromagnetic modes and their field distributions are evaluated for a dielectric 
cylindrical structure embedded in another dielectric, with a thin metallic film at the 
cylinder /dielectric interface. These modes provide energy relaxation channels for ex- 
cited dipole emitters located inside or outside the cylinder. Significantly, we find that 
the emission rate is sensitive to the magnitude of the electron density of the metallic 
film coating. For typical parameter values, we find large enhancements of the emis- 
sion rate, which can be in excess of three orders of magnitude, relative to the case in 
the absence of the film, arising at specific ranges of electron density. The theory is 
shown to conform with known limits, including the high density (perfect conductor) 
limit and the large distance (unbounded bulk) limit. The implications of the predicted 
enhancement due to the metal coating for the purpose of guiding atoms within such 
microstructures are pointed out. 

PACS Numbers: 34.50.Dy, 32.80.-t, 41.20.Jb, 42.50.-p 
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There is at present much interest in small (micrometre to nanometre) scale spher- 
ical [1,2] and cylindrical structures [3-12]. A problem of immediate interest in this 
context is the response of quantum systems to the electromagnetic fields supported by 
structures of this kind. We concentrate here on cylindrical structures, as these have 
featured more prominently in recent applications than spherical structures. Although 
some work on the coupling of fields to electric dipoles near dielectric cylinders has been 
reported recently using Green functions and other methods [8], the greater majority 
of reports have been concerned with cyhndrical structures where the walls are totally 
impenetrable to electromagnetic fields of all frequencies. In particular, Rippin and 
Knight [9] quantised the electromagnetic modes in hollow cylindrical waveguides with 
perfectly conducting walls and proceeded to calculate the decay rate of an oscillating 
electric dipole located inside the cylindrical guide. However, as we have just pointed 
out, dielectric cylindrical structures have also been the subject of investigation, espe- 
cially in experimental contexts and there is now clear need for theoretical investigations 
appropriate for the more general situation, involving dielectric structures. 

Furthermore, it has become clear that the presence of a thin metallic film as an 
overlayer on a planar dielectric structure can lead to a considerable enhancement of 
the electromagnetic fields [13-15]. For similar reasons it is envisaged that a metal- 
coated dielectric cylinder should exhibit enhanced electromagnetic fields. As far as the 
authors are aware, the effect of a thin metallic film coated at the cylinder surface on the 
electromagentic properties of such a system have not been explored before, nor has the 
coupling of the electromagnetic modes, modified by the presence of the film, to quantum 
systems localised in the vicinity of the film. In particular, the characteristics of dipole 
emission and the parameters which control the enhancement need to be understood 
and quantified. The primary aim of this communication is to report the results of a 
study along these lines. 

We consider a general system comprising a dielectric cylinder of radius Rq and 
dielectric constant ei embedded in a different material of dielectric constant 62- A 
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thin film of a given metal is introduced at the cylindrical interface between the two 
dielectrics. The system is schematically shown in Fig. 1(a) in which an emitter of electric 
dipole moment vector and oscillation frequency u>o can be situated at arbitrary points 
inside the cylinder (r^ < Rq) or outside it (rj > a) where Rq is the cylinder radius and 
rt is the cylindrical polar radial position of the dipole emitter. 

The allowed electromagnetic field modes are evaluated using the photon tunneling 
technique [16]. Briefly, the photon tunneling procedure involves plane waves which can 
either be transverse magnetic (TM) or transverse electric (TE) incident within region 2 
(outside the cylinder) at an azimuthal angle 0o- This is reflected and transmitted at the 
cylinder surface, producing appropriate linear combinations of TE and TM solutions at 
the different regions in a manner satisfying electromagnetic boundary conditions. For 
this, one needs to express the incoming plane wave form in terms of Bessel functions 
using the following identity in cylindrical coordinates r = (r^, 4>, z) [16] 

oo 
n=— oo 

where /c2 is one of two wavevectors in the problem defined by 

k, = (k,„fc,); fc,2 = ^ (^ = l,2) (2) 

In addition to the continuity of the tangential components of the electric field vector 
at the cylinder interface, the second set of boundary conditions must ensure that all 
tangential magnetic field components experience a jump due to the conductivity a 
associated with the metallic film coating 

a — — TT 

m*{uj + I'y) 

Here is the two-dimensional electron density of the film, e and m* are the electron 
charge and effective mass, respectively, and 7 << a; represents a small plasma loss 
term. This boundary condition can be expressed as [13] 

vt X [H(i?o+,0,-2) -H(i?o-,</',-2)] =^TE||(i?o,0,-2) (4) 
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where the magnetic field vectors are evaluated at -Ro± where the ± subscripts imply that 
the radial coordinate approaches rj = Rq from outside the cylinder (+) and inside it 
(— ). The resulting field distributions inside and outside the cylinder emerge in terms 
of a sum over Bessel functions, refiecting the cylindrical symmetry of the structure. 
There are four types of mode: propagating {kti and kt2 both real), evanescent {kn 
imaginary, kt2 real), guided {kn real, kt2 imaginary) and surface modes {k^ and kt2 both 
imaginary), whose characteristics depend on the relative magnitudes of the dielectric 
functions ei and 62, as well as the electron density Ug of the metallic coating. Figure 
1(b) schematically shows the spatial dependence of the four types of mode mentioned 
above in the region of the cylinder surface. Figure 2 displays the dispersion regions of 
the allowed modes in the u versus k^ plot. The light lines corresponding to ei and 62 are 
also shown. In general, the modes will have a mixture of both TE and TM characters. 
However, in the limit tt-^ — * 00, corresponding to the perfect conductor case, they split 
into pure TE and pure TM forms [16]. 

The analytical expressions of the electric field distributions corresponding to the 
allowed modes are too complicated to be displayed in full in this communication. Here 
we give a brief description of the main features with reference to the z-component of 
the electric field vector of the evanescent mode. The spatial dependence of this field 
vector component, both inside and outside the cylinder, can be displayed in the form 



£e,zir,'k2,t) = 



{kt2rt)0{rt - Rq) , 
OinJn{kt2rt)0{rt - Rq) H T~rr~\ ^ 



7nIn{i^tirt)0{Ro - n) 



^i[n{4>-<j>o-n/2)+kzZ-u)t] ^g-j 
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where 6 is the step function and J„, Hn and Kn are Bessel functions, while na = ikn 
with kti given by Eq.(|^). The subscript e in the electric field z-component in Eq.(^ 
specifies the type of mode as the evanescent mode. The forms of the factors a„, 7^ 
and An are determined by the field normalisation requirements and the boundary 
conditions, the latter giving /3„ and 7^ in terms of Note that the relation A„ = 



constitutes the dispersion relation for the surface and guided modes when a„ = and 
7„ is given in terms of /3„. Similar expressions exist for the other components of the 
field corresponding to this type and for the other types of mode. Formally, we write 
for the zth component of the quantised electric field vector 

E,{r, t) = J2 {^r,,i(r, k2, t)a(k2, v) + h.c.} (6) 

k2,»? 

where rj specifies the type of mode {rj = e and i = z m Eq.(^)); a(k2, rj) is the annihila- 
tion operator for this type of mode. In the case of propagating and evanescent modes, 
the process of normalisation for the quantised fields is carried out for a plane wave in 
an infinite bulk of material 2 (i.e. as though the cylinder does not exist) to give a„. 
For the surface and guided modes the quantisation process follows the familiar path 
to give j3n- It can be shown that the quantisation procedure amounts to determining 
normalisation factors from the condition 

eo / rf'r£,(r, k2).£;(r, ks) = \hu{k,, 77) (7) 

where £,,(r,k2) are the electric field vector functions associated with the quantised 
field mode of type ?], frequency u>{k2,ri) and wavevector k2. 

The energy relaxation rate for the oscillating dipole is evaluated by application of 
the Fermi golden rule 

r(r) = |J E {0} |-M-E(r)| {k2, v})f S{uj, - u) (8) 

where E(r, t) is the quantised electric field vector operator at the position of the electric 
dipole. The notation is such that |e) and \g) are the two quantum states representing 
the dipole system, |{0}) stands for the corresponding zero photon (vacuum) field state, 
I {k2, 7]}) is a single quantum field state of frequency a;(A;2, yf) and r] runs over the allowed 
types of mode. 

The result emerging from Eq.(|^) should be modified to take account of the local 
field correction [17] such that 



where 



e(r, u) = ei9{Ro - n) + e29{rt - Ro) 



(10) 



The expressions arising from Eq. (^ , together with Eq. (|^) , by necessity, require nu- 
merical evaluations, which permit the exploration of the changes in the characteristics 
of the system with varying cylinder radius, dipole frequency, dipole position inside and 
outside the cylinder and with varying metallic film electron density. The results are 
shown in Figs. 3 to 5. 

Figure 3 exhibits the variations in the rate with the dipole radial position for 
a cylinder of fixed radius Rq = 500nm. The inset to the figure shows the variations 
corresponding to the limiting cases = (absence of the film) [8] and Us = oo 
(perfectly conducting film) [9-11]. The values at = are for dipole position on the 
cylinder axis. On the opposite side, for rt large (r^ >> Rq), it is seen that all the 
curves converge to the same value corresponding to the rate of a dipole embedded in 
the infinite bulk of dielectric 2. In the vicinity of the film where ^ Rq marked 
variations can be seen, depending on the density of the film. For Ug = 10^°m~^ there 
is considerable enhancement near the surface; at this density the dipole frequency is 
at resonance with that of a surface mode. This is followed by a marked reduction at 
Hg = lO^^m"^. Further increase in the density to Ug = lO^^m"^ results in the rate 
dipping dramatically to smaller values approaching zero at the surface, as one would 
expect in the perfect conductor (large density) limit, shown in the inset to Fig. 3. Note 
that, for this cylinder radius, in the large limit the lowest order (n = 0) mode is the 
only mode allowed in a cylinder of impenetrable walls. 

Figure 4 displays the variations of the emission rate with the film density Ug for a 
dipole oriented parallel to the axis. The oscillation frequency is such that Hluq = 0.5eV 
and the radius of the cylinder is taken to be Rq = 500nm. Note the sudden onset of 
the dip in the rate for the dipole at the centre of the cylinder and close to the surface 
at rt = 0.9Ro. This sudden dip with increasing density is for all dipoles located inside 
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the cylinder and which oscillate at fiuo — 0.5eV. This feature can be explained by the 
fact that at such a low frequency the lowest order mode cannot be excited since the 
radius of the cylinder is below the cut-off condition Aq > Ro that exists for a cylinder 
surrounded by a perfect conductor. 

By contrast, Fig. 5 displays the variations corresponding to the case in Fig. 4, but 
for a dipole oscillating at frequency such that fkoo = 1.5eV. Note the appearance of 
the second peak at high density when the dipole is inside the cylinder near the surface 
(dashed curve) and also the finite values at high densities for a dipole inside the cylinder 
(dashed and sohd curves) . The finite values at high densities can be explained by the 
fact that the lowest order mode in the corresponding perfect conductor (cylindrical 
waveguide) case has a frequency below ojq and this provides the necessary decay channel 
leading to a finite emission rate. The second peak in the dashed curve (dipole near the 
surface) is due to the next order mode in a cylinder with impenetrable walls which, 
although it has a frequency just above cuq, is at resonance inside the cylinder with the 
incoming TM wave within a certain finite range of surface densities. 

From these results one clearly concludes that the presence of the metaUic film 
modifies the electromagnetic properties of the dielectric structure in a manner strongly 
dependent on the electron density of the metallic film. In particular, it can lead to a 
considerable enhancement of the spontaneous rate. The theory also reproduces various 
appropriate limits, most notably: the high Ug limit, corresponding to the perfect con- 
ductor case, the zero Us limit, corresponding to the bare dielectric cylinder case and 
the limit when the dipole is far away from the cylinder, corresponding to an unbounded 
medium 2. 

Cylindrical dielectric structures, are particularly interesting for use not just as elec- 
tromagnetic waveguides [3-5], but also as atom guides, where the guiding mechanism is 
mainly controlled by excited cavity modes [6,7,9,10]. It is envisaged that the develop- 
ment of atom guides at such a small scale would lead to much desirable advancements 
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in atom lithography and should facilitate atomic physics research. As mentioned previ- 
ously, the enhancement of the evanescent mode due to the introduction of the metallic 
film has been put to good use in the context of atom mirrors where an overlayer is 
deposited on the planar dielectric surface of the conventional layer structure. This has 
been shown, both experimentally and theoretically, to lead to a pronounced enhance- 
ment of the repulsive potential arising from the modified evanescent mode [12-14]. 

Similarly, the enhancement found here in the case of a metal-coated cylinder is 
envisaged to produce a repulsive force which acts on the atom near the film in a hollow 
cylindrical structure. This should lead to an efficient atom guide in such a manner 
which is dependent on the density of the metallic film. The role of the metallic film in 
aiding the guiding action in hollow metal-coated dielectric cylindrical guides is currently 
under investigation and the results will be reported in due course. 



8 



References 

1. H. T. Dung, L. Knoll and D.-G. Welsch, Phys. Rev. A 64, 013804 (2001) 

2. M. A. Kaliteevski, S. Brand and R. A. Abram, J. Mod. Optics., in press (2001) 

3. Y. Jiang and J. Hacker, Appl. Optics 33, 7431 (1994) 

4. M. A. Kaliteevski, R. A. Abram, V. V. Nikolaev and J. S. Sokolovski, J. Mod. 
Optics, 43, 875 (1999) 

5. A. J. Ward and J. B. Pendry, J. Mod. Opt. 43, 773 (1996) 

6. J P Bowling and J Gea-Banacloche, Adv. Atom Mol. Opt. Phys. 37, 1 (1997) 

7. V. I. Balykin, Adv. Atom. Molec. Opt. Phys. 4, 181 (1999) 

8. N. Nha and W. Jhe, Phys. Rev. A 56, 2213 (1997); W. Zakowicz and M. 
Janowicz, Phys. Rev. A 62, 013820 (2000) 

9. M. A. Rippin and P. L. Knight, J. Mod. Opt. 47, 807 (1996) 

10. S. Al-Awfi and M. Babiker, Phys. Rev. A 58, 2274 (1998) 

11. M. Babiker and S. Al-Awfi, J. Mod. Optics 48, 847 (2001) 

12. D. Miiller, E. A. Cornell, D. Z. Anderson and E. R. I. Abraham, Phys. Rev. A 
61, 03341 (2000). 

13. C R Bennett, J B Kirk and M Babiker, Phys. Rev. A 63, 033405 (2001) 

14. S Feron, J Reinhardt, S Lebouteux, O Gocreix, J Baudon, M Ducloy, J Robert, 
C Miniatura, S N Chormaic, H Haberland and V Lorent, Opt. Commun. 102, 

83 (1993) 

15. T Esslinger, M Weidenmiiller, A Hammerich and T W Hanch, Opt. Lett. 18, 
450 (1993) 

9 



16. A. Ishimaru, "Electromagnetic wave propagation, radiation and scattering" (Pren- 
tice Hall: New Jersey, USA 1991) 

17. R. J. Glauber and M. Lewenstein, Phys. Rev. A 43, 467 (1991) 



10 



Figure Captions 



Figure 1 

(a) A schematic drawing showing the dielectric cyhnder of material 1 immersed in 
another dielectric material 2. A thin metallic film is coated at the cylinder surface. 
The dielectric constants are ei and 62 and the dipole is denoted by an arrow, (b) 
Schematic representation of the spatial dependence of the four types of mode (propa- 
gating, evanescent, guided and surface modes) in the region of the cylinder surface. 

Figure 2 

The regions of dispersion of the allowed modes in the cyhnder. The light lines for the 
bulk media 1 and 2 are shown for the ei = 1 and €2 = 4 case. The left shaded region 
contains the propagating modes, the other shaded region contains the evanescent modes 
and the solid curves are the lowest orders of an infinite number of surface modes, the 
lowest in energy being the n — mode and the arrow signifying increasing n. 

Figure 3 

Variations in the spontanous emission rate with the dipole radial position for a 
cylinder of fixed radius Rq = 500nm. The oscillation frequency is taken to be such 
that hujo — 1.5eV and the dipole moment vector is oriented parallel to the cylinder 
axis. The other parameters are €1 — 1 and €2 = 4. The three curves correspond 
to different metallic film coating density: = 10^°m~^ (solid curve), rig = lO^^m"^ 
(dashed curve) and Ug — lO^^m"^ (dot-dashed curve). The inset to the figure shows 
the variations corresponding to the limiting cases = 0, absence of the film (solid 
curve), and Us = 00, perfectly conducting film (dashed curve). 
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Figure 4 



Variations of the emission rate with the film density Us for a dipole oriented parallel 
to the axis. The oscillation frequency is such that Hcuq — 0.5eV and the radius of the 
cyhnder is taken to be Ro — 500nm. The three curves are for three different positions 
of the dipole: = 0, dipole at the centre (solid curve); = 0.9i?o, dipole inside, 
near the film (dashed curve); and = 1.5i?0) dipole outside the cylinder (dot-dashed 
curve) . 

Figure 5 

Variations corresponding to the case in Fig. 4, but for a dipole oscillating at frequency 
hcuo — 1.5eV. The three curves correspond to three different positions of the dipole: 
rt — 0, dipole at the centre (solid curve); rt — 0.9Ro, dipole inside, near the film 
(dashed curve); and rt — 1.5Ro, dipole outside the cyhnder (dot-dashed curve). 
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Figure 1a Banyard 
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Figure 1b Banyard 



14 



15 




Figure 3 Banyard 
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Figure 4 Banyard 
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Figure 5 Banyard 
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